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Abstract 

It is shown that three SU(2) singlet neutral scalars (two CP-even and one CP-odd) in the 
spectrum of models based on the gauge symmetry SU (3) c <g> SU (3)/, (g) U(l)x , which do not contain 
exotic electric charges, are realistic candidates for thermally generated self-interacting dark matter 
in the Universe, a type of dark matter that has been recently proposed in order to overcome some 
difficulties of collisionless cold dark matter models at the galactic scale. These candidates arise 
without introducing a new mass scale in the model and / or without the need for a discrete symmetry 
to stabilize them, but at the expense of tuning several combinations of parameters of the scalar 
potential. 

PACS numbers: 95.35.+d, 12.60.Fr, 14.80. Cp 
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The inflationary model plus the cold dark matter model supplemented by a cosmological 
constant have been very successful in explaining the observed structure of the Universe on 
large scales. Notwithstanding, in the past years very high resolution N-body simulations have 
shown that collisionless cold dark matter (CDM) models produce far too much substructures 
on galactic scales as compared with observations and predict dark matter halo density 
profiles that diverge at the centre in disagreement with the constant density cores observed 
in late- type dwarf and low surface brightness galaxies j^. 

A possible way out to these problems is to assume that the dark matter particles are 
self-interacting with a large scattering cross-section but negligible annhilation or dissipation 
P). However, self-interacting CDM models lead to spherical halo centers in clusters which 
is not in agreement with ellipsoidal centers indicated by strong gravitational lensing data 
[J] and have difficulties to reconcile the number of satellites in the Galactic halo with the 
observed number of dwarf galaxies in the Local Group {5j. In spite of these difficulties, 
self-interacting dark matter models are well-motivated as alternative models. 

The unsatisfactory ad hoc addition of just one singlet neutral scalar field, and a discrete 
symmetry in order to guarantee its stability, as the simplest possible renormalizable exten- 
sion of the Standard Model that includes non-baryonic self- interacting dark matter jfJQ], has 
motivated several authors to investigate the possibility for such a class of dark matter arising 
naturally from realistic gauge models. In particular, Refs. and |9( have studied this issue 
in the context of the Standard Model (SM) extension based on the SU(3) c ®SU(3)l<£)U(1)x 
gauge symmetry (the 3-3-1 model). In Ref. js], by introducing a new very small energy scale 



10 GeV), a physical CP-odd neutral field satisfaying the Spergel-Steinhardt bound 



3] 



for self-interacting dark matter has been identified in the scalar sector of the Pisano-Pleitez- 
Frampton model (PPF) Q], which is a model that includes quarks with exotic electric 
charges —4/3 and 5/3 and double charged Higgs and gauge bososns. In Ref. [I| and after 
imposing by hand a discrete symmetry, two scalar fields, one CP-odd and one CP-even, have 
been identified as realistic candidates for self-interacting dark matter in the version of the 
model referred to as "3-3-1 model with right-handed neutrinos" which is a model where 
the particles have only ordinary electric charges. 

Our goal in this letter is to search for realistic candidates for thermally generated self- 
interacting dark matter arising from the scalar sector of 3-3-1 models without exotic electric 
charges, without the need of imposing an extra symmetry and/or without introducing a 
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new mass scale, and taking care of the consistency of the full scalar spectrum of the theory. 
Thermally generated singlet scalars with MeV range masses and 0{Mpl) couplings to matter 
was proposed for the first time as self- interacting dark matter in p, and has been recently 
proposed as a candidate for explaining the observed fluxes of 511 KeV photons from the 



Galactic bulge [12j . Clearly, this range of masses is not in the 1 GeV - 1 TeV range preferred 
by most models of weakly interacting massive particles (WIMPS) 

The 3-3-1 extension of the SM has received special attention over the last decade since 
several models can be constructed so that when anomaly cancellation takes place between 
the fermion families (three-family models), and not family by family as in the SM (one- 
family models), these models predicts that the number Nj of families must be an integer 
multiple of the number N c of colors of SU(3) C . Moreover, since SU(3) C asymptotic freedom 
requires N c < 5, it follows that Nf = N c = 3. In this way these models provide a possible 
answer to the question of the number of fermion families in nature 1G, 141 • 

All possible 3-3-1 models without exotic electric charges are presented in Ref. |15l |. These 
models share in common not only the same gauge boson sector but also the same scalar sec- 
tor. For the 3-3-1 models both, SU(3)l and U(l)x are anomalous, so special combination 
of multiplets must be used in each particular model in order to cancel the several possible 
anomalies, and end with physical acceptable models. It is shown in Ref. that by demand- 
ing cancellation of the triangle anomalies: [SU(3) L ] 3 , [SU(3) C } 2 U(1) X , [SU(3) L } 2 U(1) X , 
[grav] 2 U(l)x (the gravitational anomaly), and [£/(l)x] 3 , only ten different models with- 



out exotic electric charges can be constructed. Two of them are one-family structures and 
have been studied in Refs. and ^J. The other eight models are three-family models. 
One of them is the already mentioned "3-3-1 model with right-handed neutrinos" [ill, an- 

n n 

other one has been analyzed in |18j . two more have been studied in [19(, and the remainig 
four models have not yet been studied in the literature, as far as we know. 
Now, if our aim is to break the symmetry in the way 

SU{3) e ®SU(3) L ®U(l) x 

— > SU($) c ®SU(2) l ®U(1)y 
— > SU(3) C ®U(1) Q , 

(1) 

and at the same time give masses to the fermion fields in the model, the following three Higgs 



scalars must be introduced [llj]: $i(l, 3*, — 1/3) = , 05,0']°) with Vacuum Expectation 
Value (VEV) (^} T = (0,0, V); $ 2 (1, 3*, -1/3) = (02,02,02°) with VEV ($ 2 > T = (0,^,0), 
and $ 3 (1, 3*, 2/3) = (0^, 0^ , 3 + ) with VEV ($ 3 ) T = (w,0,0), with the hierarchy V > 
v ~ u ~ 174 GeV, the electroweak breaking scale. The consistency with the low energy 



phenomenology requires \/u 2 + 1> 2 = 174 GeV (U|, and by using experimental results from 
the CERN LEP, SLAC Linear Collider and atomic parity violation data, a lower bound on 
the mass scale V have been calculated in Refs. jig, [l7| ; generically, V > 1 TeV. 

The most general renormalizable scalar potential which includes $1, $2 and $3 can be 
written in the following form 



V($i,$ 2 ,$ 3 ) = + A«2$2$2 + /is$3$3 + y($l$2 + /l.C.) 

+A 1 ($ t 1 $ 1 ) 2 + A 2 ($ 2 $ 2 ) 2 + A 3 ($ 3 $ 3 ) 2 + y [($l$ 2 )(<f t 1 $ 2 ) + /i.e.] 
+A 5 ($1$ i)($ t 2 $ 2 ) + A 6 ($ t 1 $ 1 )($ 3 $ 3 ) + A 7 ($ t 2 $ 2 )($ 3 $ 3 ) 
+A 8 ($ t 1 $ 2 )($ t 2 $ 1 ) + A 9 ($ t 1 $ 3 )($ 3 $i) + Aio($^ 3 )(^$ 2 ) 
+i[A 11 ($ t 1 $ 1 )($ t 1 $ 2 ) + A 12 (^$ 2 )($ t 1 $ 2 ) + A 13 ($ 3 <l>3)($ t 1 <l> 2 ) 
+Ai 4 ($ t i$3)($ 3 < f2) - fe ijk <S>\& 2 <$> k 3 + h.c], (2) 



which is the same scalar potential studied in Refs. [20j, |2l|. In this equation the fi's and / 
are coupling constants with dimension of mass and the A's are dimensionless. 

For the sake of simplicity we assume that the VEV are real and, for convenience in 
reading, we rewrite the neutral scalar fields as 



,0 <PiR + i<Pii 



'1 



V2 



= v+ ^ + 



V2 



L0 _ „ , 02fl + #2/ 

j 2 — y -r ■ 



M + rk±^3/ , (3) 



In the literature, a real part 0r is called a CP-even scalar, and an imaginary part 0/ is 
called a CP-odd scalar or pseudoscalar field. 
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A complete study of the scalar mass spectrum for the type of 3-3-1 models we are inter- 
ested in has been done in Ref. jsjj, in the limit / ~ V ^> u,v. In this limit, the couplings 
of one of the 3-3-1 CP-even scalars to the SM fermions and gauge bosons become identical 
to the couplings of the SM Higgs, and the mass spectrum for the scalar fields both in the 
charged and in the neutral (CP-even and CP-odd) sectors, is as follows (2]]. 

For the charged scalar sector the square mass matrix has two eigenvalues equal to 
zero equivalent to four would be Goldstone bosons (Gf,Gf ). Now, with C a = cos a = 
uj \Ju 2 + v 2 , S a = sin a = v/y/u 2 + v 2 , and in the limit / ~^>u,i;, the physical charged 
scalars are 

Hf ~ C a( pf + S a <f)i, (4) 

with squared masses at the TeV scale given by 

m 2 H± ~ \ 9 V 2 - lfV V -, 
w i 2 u 

m% ± * \fV (- + -), (5) 

respectively. 

For the CP-even scalars the square mass matrix has one eigenvalue equal to zero corre- 
sponding to another would be Goldstone boson G\ ~ —<Pir- Defining 

2(A 8 + A 4 -4Ai)i; + m 



cos 2/3 
sin 2(3 



[2(A 8 + A 4 - 4Ai)« + u} 2 + I6\ 2 u v 2 
AX n v 



'[2(A 8 + A 4 - 4Ai)v + u] 2 + I6\ 2 u v 2 
and in the limit / ~ V ^> u, v, the physical CP-even scalars and their squared masses are 

Hi - CA - Satin rn 2 H o ~ ±fV; 

Hi ~ S a (f 2R + C a (f> 3R , m 2 H o ~ uv; 

H° 3 * C^ R - S^' 2 ° R , m 2 Hi ~ ^V 2 ; 

Hi ~ 5^ + m|o ^ ^A 2 ^ 2 , (6) 

where 



A a , 2 = A + ± + 4A?! > 0, 

A ± = (/n/2^) + A 4 + A 8 ±4A 1 , (7) 



and Cp, Sp stand for cos/3, sin/3, respectively. 

Notice that H$ and H% are mainly mixtures of SU(2) singlets and that their masses 
depend on the combinations of parameters Ai and A 2 which are not fixed by the experiment, 
except for \/u 2 + v 2 = 174 GeV. So, m H o and m H o can in principle be tuned in order to 
have two light CP-even scalars in the low energy spectrum of the model. Notice also that 
H® and H 2 are mainly mixtures of members of SU{2) doublets. Besides, it has been shown 



in Ref. |2l| that # 2 ° is to be identified with the SM Higgs [See Eq. fljft below]. 

For the CP-odd scalars the square mass matrix has three zero eigenvalues associated 
with the three additional would-be Goldstone bosons: G\ ~ G® ~ 0^, and G® ~ 
— S a (f>2i + Cq03/- The nonzero eigenvalues correspond to the two physical states 

h\ ~ c a <p° 2I + s^, 

hi ~ (8) 



with squared masses given by 



1 ,u 2 + v 2 



m h° ~ if v > 

ftl 4 uv 



1 [fuV - 2v(X 4 - \ 8 )V 2 } 

(9) 



4 v 
1 [u - 2v(X 4 - X 8 )]V 2 
4 v 



respectively. 

Also in this sector, the scalar h® is mainly a SM singlet and its mass depends on the 



difference u — 2w(A 4 — A 8 ) , where yvF+v 2 = 174 GeV. Since the value of A4 — As is not fixed 
by the experiment, we can tune the value of m h o so that we also have a very light CP-odd 
scalar in the low energy spectrum of the model. Notice that h\ is heavy and is mainly a 
mixture of electrically neutral components of SU (2) doublets. 

It is worth to say that, as required by the consistency of the model, the original eighteen 
degrees of freedom in the scalar sector have become eight would be Goldstone bosons (four 
electrically neutral and four charged), and ten physical scalar fields, six neutrals and four 
charged ones. 

Since the CP-even scalars H® and H% and the CP-odd scalar h® are mainly SM singlets, 
they couple very weakly to the scalars h\, Hq and H® (which are mixtures of electrically 
neutral components of SU (2) doublets) and should couple very weakly to the SM fermions 



and gauge bosons. In fact, the trilinear couplings between the physical scalars of the model 



and the SM gauge bosons have been calculated in Ref. 
in the limit / ~ V ^> u, v , is 



2l| and the corresponding Lagrangian, 



Linear = (q M W± W»~W+ + ^-M Z Z» Z H° 2 

-[(eM w± # + gMzS&Z^W+Gl + h.c 
+i(p - rfW+H^hl + h.c.V 



(10) 



where the scalars h®, and H® are absent and from which it is clear that the couplings of 
H® coincide with the ones of the SM Higgs. 

These features of h®, and H® make them stable without the need of imposing an extra 
symmetry, and make also them suitable candidates for self-interacting dark matter in the 
Universe. Let us see. 

A realistic candidate for self-interacting dark matter must have mean free path A, at 
the solar radius in a typical galaxy, in the range 1 kpc < A < 1 Mpc, where A = l/(na). 
n = p/m^a is the number density of a generic singlet scalar field h°, a is the cross section 
for the self-interaction h°h° — > h°h°, and p is the mean density of a typical galaxy (p = 0.4 
GeV/cm 3 jsj]). The bound on A can be translated into a bound on the quotient a/m h o, 
namely 

2 x 10 3 GeV~ 3 < — < 3 x 10 4 GeV~ 3 . (11) 

Near the threshold (s ~ 4m^ ) the cross section o is roughly given by a = A|/(167rs) ~ 
A|/(647rm^ ), where As is the coupling constant for self-interaction. The value of As for h®, 
H® and H®, obtained from the scalar potential written in terms of physical fields, are given 
respectively by 

Asi = A 2 , 

^S2 = Ai cos 4 (3 + A 2 sin 4 (3 

+^(A 4 + A 5 + \ 8 )sin 2 2p 

— -(An cos 2 (3 + A12 sin 2 0) sin 2/5, 
— Ai sin 4 (3 + A 2 cos (3 
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+-(A 4 + A 5 + \ s )sin 2 2f3 

+— (An sin 2 (3 + A12 cos 2 (3)sin2(3. 



(12) 



Let us write the mass m^o of each one of our dark matter candidates in the form 



rrii = 




(13) 



where i = 1,2,3, A 1 = \u - 2v(X 4 - A 8 )|, A 2 = v|Ai|, and A 3 = v|A 2 |, for h° 2 , H$ and Hi, 
respectively, with A ij2 as defined in Eq. (J7J). 

With Xsi = 0.5 (a natural value as compared to the SM Higgs self-coupling), V = 1 
TeV, v = 123 GeV, and if we tune each one of the Aj so as to lie in the range 5.90 x 10 -9 
GeV < Aj < 3.57 x 10~ 8 GeV, we obtain the required Spergel-Steinhardt bound in Eq. (fTTj) 
and, for each one of the dark matter candidates provided by the present model, a tree-level 
mass in the range: 3.46 MeV < 8.52 MeV, a typical value for thermally generated 

self- interacting dark matter scalars 6]. The range of values for each Aj can be tuned with 
a variety of values of the parameters involved in their definitions. 

The scalars h®, and with a mass in the range of a few MeV, are non-relativistic 
in the decoupling era (T dec ~ 1 eV), then, for a SM Higgs boson with a mass of the order of 
100 GeV the dominant contribution to the cosmic density of these scalars comes from the 
decay of thermal equilibrium scalar if 2 (the scalar of the model to be identified with the 
SM Higgs) to h°h pairs at temperatures less than the temperature Tew of the electroweak 
phase transition (Tew > l-5m#o) 6]. The trilinear couplings of /i^, -£^3 and to if 2 , also 
obtained from the scalar potential in Eq. (0) written in terms of physical fields, have the 
following strengths: 
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+ ^ cos a{2A 6 sin 2 (3 + 2A 7 cos 2 (3 

+A 13 sin2/3}. (14) 

The cosmic density of light gauge singlet scalars from the H® decay has been calculated in 
Ref. 0] and is given by 

n h o = 2g(T,)T 3 ^^, (15) 



with 



0. = * = V T i (16) 
l "T3 An 3 Km 3 ' 1 j 

where T 7 = 2.4 x 10~ 4 eV is the present photon temperature, g(T 7 ) = 2 is the photon 

degree of freedom, g(T) = gs + 7gi?/8 (gs and gp are the relativistic boson and fermion 

degrees of freedom, respectively), p c = 7.5 x 10~ 47 h 2 GeV 4 is the critical density of the 

Universe (h ~ 0.71 is the Hubble constant in units of 100 km s _1 Mpc -1 ), rj = 1.87, 

K 2 = 47r 3 #(T)/45m 2 / and m pl = 1.2 x 10 19 GeV is the Planck mass. We will take T = m H o, 

since most of the contribution to each 0j comes from T < m H o < Tew Q|- 

The value of g(m H o) depends on the particle content of theparticular 3-3-1 model under 

study, but typically g{m H o) ~ 130 for a three-family model jq. For simplicity we take all 

the rrii equal to a common value m = 5 MeV, and m H o = 150 GeV. Then, from Eqs. ()15j) 

and (fTHjl we get Q h o ~ 0.3 for J2i T 2 — 1-88 x 10 -12 GeV 2 . Assuming that each Tj contributes 

with approximately the same value to X^r 2 , we obtain Tj ~ 10~ 6 GeV, in broad agreement 



with the estimates in 



implying that such singlet scalars are very weakly coupled to the 
SM sector. Also in this case, as can be seen from Eq. (|14jh this value for J2i T 2 can be tuned 
with a variety of values of u and v and of the A's involved, a task that can be done more easily 
by taking, for example, A 7 < without spoiling the scalar mass spectrum. Consequently, 
the three light singlet scalars h®, and H® do not overdose the Universe. 

In conclusion, for 3-3-1 models without exotic electric charges we have considered the 
scalar sector and we have shown that for all these models this sector provides three light 
SU(2) singlet neutral Higgs bosons, two CP-even and one CP-odd, each one with a mass in 
the range of a few MeV, and each one satisfying the properties of a candidate for thermally 
generated self-interacting dark matter in the Universe, all of them at the expense of tuning 
the values of combinations of several of the parameters in the scalar potential. These scalars 
have feeble interactions with ordinary matter, do not overpopulate the Universe and are 
stable without the need of an extra discrete symmetry. 
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The same problem was studied in Ref. |8| for the PPF model, where the authors identify 
a light CP-odd scalar as a dark matter candidate at the price of introducing a new very 
small energy scale, different from V, u and v, and of the order of 10~ 7 GeV, associated to the 
trilinear coupling / in the scalar potential. A related analysis is done in |9( for the so-called 
"3-3-1 model with right-handed neutrinos" (a model without exotic electric charges) when 
the discrete symmetry $1 — > — $1 is imposed, with the result that the model provides two 
scalar dark matter candidates, one CP-even and one CP-odd, also at the price of fine-tuning 
the parameters of the scalar potential. However, it is clear from Eq. (J2J) that this discrete 
symmetry forbids, among others, the trilinear coupling /, which plays a fundamental role 
in our analysis. So, what we have done is to generalize the study presented in Ref. [I| by 
considering the most general scalar potential produced by the three Higgs scalar triplets of 3- 
3-1 models without exotic electric charges, without the introduction of an extra symmetry to 
account for the stability of the dark matter and that drastically reduces the scalar potential, 
and avoiding the introduction of a new mass scale. We also want to emphasize that our 
analysis is valid for all the 3-3-1 models without exotic electric charges, and not only for the 
so-called "model with right-handed neutrinos" . 

The main conclusion of this study is that models based on the 3-3-1 gauge structure can 
provide good candidates for thermally generated self-interacting scalar dark matter but at 
the price of: introducing a new mass scale as done in Q, or introducing an ad- hoc discrete 



a, 



symmetry as done in or by tuning some parameters of the scalar potential as shown in 
this letter. 
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